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surface displacements resulting from the detonation of a small 
explosive charge placed on the object surface. The primary wave- 
fronts resulting from the explosive charge were sufficiently 
large as to be visible upon reconstruction of the holographic 
interferograms. In addition, details of the surface wavefronts 
which could not be deduced using conventional holographic inter¬ 
ferometry were measured and displayed by modifying the pulsed 
recording geometry and subsequently applying heterodyne inter¬ 
ferometric analysis to the holograms to interpolate between the 
observed interferometric fringes. Although theory predicts that 
still greater resolution is possible, fringe interpolation to 
1/900 of a fringe (about 3 a displacement) was demonstrated. 

Based on these experimental results and the mathematical rela¬ 
tionships which have been developed to describe system perfor¬ 
mance, it appears likely that a practical system for studying 
sub-Angstrom transient phenomena may be realized using 
heterodyne holographic interferometry. 












Abstract 


A technique for recording and visualizing acoustic wavefronts over a finite field 
of view and with a resolution on the order of Angstroms has been developed. 
Multiple exposure holograms of an aluminum test block were used to record the 
surface displacements resulting from the detonation of a small explosive charge 
placed on the object surface. The primary wavefronts resulting from the explosive 
charge were sufficiently large as to be visible upon reconstruction of the holograph¬ 
ic interferograms. Tn addition, details of the surface wavefronts which could not 
be deduced using conventional holographic interferometry were measured and 
displayed by modifying the pulsed recording geometry and subsequently applying 
heterodyne interferometric analysis to the holograms to interpolate between the 
observed interferometric fringes. Although theory predicts that still greater resolu¬ 
tion is possible, fringe interpolation to 1/900 of a fringe (about 3 a displacement) 
was demonstrated. Rased on these experimental results and the mathematical rela¬ 
tionships which have been developed to describe system performance, it appears 
likely that a practical system for studying sub-Angstrom transient phenomena may 
be realized using heterodyne holographic interferometry. 
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INTRODUCTION 


To date, conventional techniques for measuring and studying transient surface 
wave phenonena have been limited to transducer techniques which provide a time 
record of surface displacement at a single position on an object surface. Piezoelec¬ 
tric, capacitive, magnetic, and single point interferometric techniques all fall into 
this category. For many testing and measurement applications, these techniques 
are quite adequate. However, some studies are difficult if not impossible to per¬ 
form with such transducer systems. If one wishes to investigate the interaction of 
surface waves with a surface defect, for example, a single point transducer must be 
moved around the defect area and the experiment repeated for each new transduc¬ 
er location in order to map out changes in the amplitude and shape of the surface 
wavefront as it is diffracted and attenuated by the defect. For some transducer 
types the effects of mounting the transducer on the surface must also be taken into 
account Further, the acoustic wave source must be very reproducible so that one 
is assured that for each new placement of the transducer, the wave which is pro¬ 
pagated is identical to the one before it Alternatively, one can imagine an experi- 
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ment where an array of noncontact ing transducers are positioned on an object sur¬ 
face so that surface displacement over a fixed field may be recorded as a function 
of time and position for a single acoustic event. The cost of implementing such a 
system and restrictions on object surface geometry would seem to preclude the use 
of such a system as a practical measurement technique. 

Holographic techniques, on the other hand, provide a means to interferometri- 
cally compare changes in surface topography of objects of arbitrary shape over a 
finite field of view. Tune-averaged holographic interferometry has been used to 
study the rather large surface displacements associated with mechanical resonances. 
Transient surface displacements associated: with impact testing and high energy ul¬ 
trasonic testing have been recorded using double exposure pulsed holography. In 
both of these holographic techniques, inierferometric fringes are observed in the 
reconstructed holographic image wherever the surface is displaced by an amount 
greater than one half the wavelength of the illuminating laser light. Thus displace¬ 
ments less than about 250nm (2500A) do not give rise to holographic fringes. 
Consequently sub-fringe wavefronts or sub-fringe detail of large wavefronts cannot 
be discerned by either of these holographic techniques. However, by modifying 
the holographic recording geometry and applying heterodyne analysis techniques to 
the recorded hologram, full field measurement and visualization of sub-fringe sur¬ 
face displacements such as those associated with traveling surface acoustic wave 
may be achieved. 


Full field recording is achieved using pulsed holography to "freeze* surface mo¬ 
tion at a particular instant in time. Using double exposure holography, the surface 






- 3 - 


topography of a test specimen at different instants in time may be interferocnetri- 
csily compared. The double exposure hologram recoids two virtual images of the 
test specimen which, when reconstructed by reillumination of the hologram with 
its reference beam, appear in space at the original object location relative to the 
hologram plate. These two virtual images may be used as virtual objects for an in¬ 
terferometer if reimaged on an observation plane. The interference pattern pro¬ 
duced at the observation plane is identical to that which would be observed if the 
light reflecting from the surfaces of two nearly identical real objects were combined 
by substituting the object surfaces for the mirrors of a Michelson interferometer. 

This Michelson interferometer analogy for holographic interferometry proves 
to be useful in predicting the effects of surface displacements on the observed 
holographic interference pattern. In fact the results obtained using this analogy are 
exact for the case where both the viewing angle and the angle of object illumina¬ 
tion are normal to the object surface when the hologram is exposed and where 
plane wave illumination is used. Assuming then that this is indeed the case one 
can see that local changes in surface displacement alter the optical pathlength from 
the object surface to the point of observation. If at a single point on the object 
there is no displacement between the first and second holographic exposures, then 
the light reflecting from corresponding points of the two virtual objects travels ex¬ 
actly the same distance to the observation plane. Consequently, the electric fields 
of the two reflected rays will be in phase and interfere constructively to give a 
bright point in the observation plane. If, on the other hand, at a different point on 
the object, some displacement occurred between exposures, then the light 
reflecting from corresponding points on the two virtual objects will not be in phase 










at the observation plane. If the displacement which occurred between exposures 
was sufficient to cause an optical pathlength change of one half the wavelength of 
the illuminating light, then the two rays from the corresponding paints on the vir¬ 
tual object will differ in phase by 180* thus interferring destructively giving a dark 
point in the observation plane. As in the case of a Michelson interferometer, a- 
object motion of one-half wavelength yields one whole wavelength of optical path- 
length change. Therefore the interferometric fringes formed over the image of the 
object in the observation plane appear to map out contours of constant displace¬ 
ment where the contour interval is one half the wavelength of the illuminating 
light. 

With the establishment of this half wavelength contour, or fringe, interval the 
sensitivity limit of visually interpreted holographic interferograms is set. In the ex¬ 
periment reported in this work, an argon ion laser with an output wavelength of 
514.5nm was used. Thus one full fringe corresponds to a surface displacement of 
257.3nm. One can readily discern one-half a fringe and therefore may resolve 
down to 128.6nm by simply counting fringes observed in the output plane of the 
holographic interferometer. Unfortunately surface displacement amplitudes gen¬ 
erated by conventional ultrasonic transducers may range from 0.01 nm to tens of 
nanometers and are, therefore, well below the limit of resolution for visually inter¬ 
preted holographic interferometry. In other words, if a flat mirror wee used as a 
test object and a hologram was exposed before and during excitation of the mir¬ 
rored object by a conventional ultrasonic transducer, no fringes would be observed 
in the observation plane of the holographic interferometer. 








Sub-fringe displacement information, however, is recorded by the hologram 
and may be extracted using heterodyne interferometric techniques to read out the 
holographic record- Heterodyne techniques employ any of several methods to en¬ 
sure that light arriving at the output plane from one virtual object, is at a slightly 
different frequency from the light arriving from the other virtual object. The 
effect of frequency shifting one of the light beams is to cause the intensity of the 
image in the output plane to vary periodically at the difference frequency of the 
two beams. While the frequency of the intensity variation is the same at every 
point on the object surface, different points in the image plane may not fluctuate 
precisely in phase with each other. The degree of phase difference in intensity 
fluctuation between two poinis in the observation plane is directly related to the 
difference in surface displacement between the two points. Thus if a one-half 
fringe difference were observed between two poinis using conventional holographic 
interferometry, there would be a 180° phase difference in intensity fluctuation 
measured between the same two points in the heterodyne case. Since one fringe 
corresponds to a 360“ phase difference, then each degree of phase difference 
corresponds to an object displacement of 257.3nm x 1/360 or about 0.7nm for the 
514.5nm line of the argon ion laser. Note that sensitivities in this range should 
permit detection of displacements from conventional ultrasonic sources. The abili¬ 
ty to resolve much smaller displacements is limited by factors which may be op¬ 
timized for a particular application. These factors include quantum noise at the 
detectors, speckle noise in the image, and environmental effects. 
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Conaderation must also be given to the process of holographically recording a 
traveling acoustic wavefront. For the studies reported in this work, a pulsed 
NdiYAG laser with a frequency doubled output at 532nm was used. The laser had 
a fixed pulse duration of 9ns and a peak output energy per pulse of about 600mJ. 
The diffraction efficiency of holograms recorded of objects in motion is a function 
of the velocity of the displacements along the direction of the bisector of the il¬ 
lumination and viewing angles. For periodic displacements, this velocity is deter¬ 
mined by the product of the frequency and amplitude of the surface waves. For 
the amplitude and frequency spectrum of the surface waves studied during this in¬ 
vestigation, the 9ns laser pulse duration was sufficiently short to ensure only a 
0.5% reduction in holographic diffraction efficiency. Thus holographic reconstruc¬ 
tion efficiency is not significantly degraded by small acoustic perturbations of the 
object surface. Likewise, the accuracy with which such small perturbations are 
recorded was shown to remain sufficiently high to permit accurate surface displace¬ 


ment measurements to the sub-Angstrom level. 






BACKGROUND 


The advantages of using pulsed holography and heterodyne analysis for study¬ 
ing traveling surface acoustic waves become apparent once the operating principles 
are well understood. To obtain such an understanding it is useful to review the 
principles of classical interferometry and holography. This section includes a dis¬ 
cussion of interferometry as a tool for measuring dynamic surface displacement. 
The Michel son interferometer is used as an example with special attention paid to 
noise sources which limit displacement sensitivity and are to some degree common 
to the heterodyne holographic system. These noise sources include pathlength in¬ 
stabilities, detector shot noise, and thermal noise in amplifiers. Heterodyne inter¬ 
ferometry, the nature of heterodyne signals, and its effect on overall signal to 
noise ratio is discussed. 

Holography relieves the restrictions imposed by classical interferometry on the 
types of test objects that may be analyzed. Typically one may inspect flatness of 
mirrors or symmetry of convex or concave optical surfaces using classical inter¬ 
ferometers. Using holography, however, one may inspect more complex and even 
diffuse surfaces. A hologram may be used to reconstruct any arbitrary reference 
wavefront especially that reflected by the test object at some earlier point in time. 
This reference wavefront may then be combined with the object wavefront to pro¬ 
duce an interference pattern over the image of the object indicative of surface dis¬ 
placements which may have occurred subsequent to the initial holographic expo- 
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sure. Thus, the entire object within the field of view of the holographic inter¬ 
ferometer may be holographically recorded and interferaraetrically compared. The 
principles of holographic wavefront reconstruction, the use of pulsed recording 
techniques for freezing rapidly occurring displacements, and the application of 
heterodyne techniques to improve holographic sensitivity will be discussed. 


Interferometry 

The Michelson interferometer shown schematically in Figure 1 will be used to il¬ 
lustrate the interferometric principles applied to surface displacement measure¬ 
ments. Tn every case it will be assumed that coherent, single mode, polarized laser 
illumination is used. The laser beam entering the interferometer from the left is 
split by partial reflectance into two beams at the beamsplitter. The portion of the 
beam deflected upwards in the figure is reflected by the reference mirror back 
through the beamsplitter and onto a photodetector. That portion of the laser beam 
not initially deflected by the beamsplitter is reflected from the test object surface 
and returned to the splitter where it is combined with the reference beam and 
deflected onto the photodetector. 

The object and reference waves arriving at the detector may be expressed by 
their respective analytic signal representations (Yariv, 1976) which is the harmonic 
solution to the free space optical wave equation for the electric field at the detector 
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Figure 1: Michelson Interferometer 
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location: 


U R -A R e>« (1) 

U 0 -Ao* 


where 

U R — reference beam 
U Q — object beam 

^ — field amplitude of reference beam 
A 0 — field amplitude of object beam 
at — radian frequency of the laser 
t — time 

k — propagation constant 2 tt/A 
dl — total pathlength difference 
between the reference 
and object beams. 


These expressions represent the amplitudes of tte two optical beams. Photo¬ 
detectors, however are square-law devices sensitive to radiant intensity not ampli¬ 
tude. Therefore the electrical signal generated by the detector will be proportional 
to the light intensity incident on the detector. 


Id-\U 0 +V r I 2 
- 0 U 0 +U R Wo+U R ) 


( 2 ) 
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where 

u o — complex conjugate of U Q 
U R “ complex conjugate of U R . 


Thus, 


iD-At+Ak+AoAdet'+e-* 1 ] 
/lucos kt11 


ami 


If A o“ A R’ then 

f D -L4Hl+et» kdl). (3) 

Since in both legs of the interferometer, the beam is folded back on itself, the 
pathlength difference^, resulting from a normal surface displacement, D, is twice 
the displacement amount. Equation 3 may now be rewritten: 

Id—IAKX+cos 2 kD) 

-lAHl+tosi^D)) (4) 

It can be seen then that minima in the light intensity falling on the detector occur 
when/)— n^-, n — 1,2,3... 


Although the Michelson interferometer in this basic form has been used to 
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detect acoustic signals (Wellman, 1980), several modifications may be made to im¬ 
prove its performance. By twice differentiating Equation 4 with respect to D one 
observes that the greatest change in intensity for a given surface deflection less 
than A/2 occurs when D — (2n+l)A/8. One approach, then, to obtaining max¬ 
imum sensitivity to surface displacements is to actively stabilize the optical path 
difference between the object and reference legs using an electro-mechanical feed¬ 
back arrangement to drive the reference mirror. An additional benefit realized by 
using an elec tronmec ham cal feedback system is that the surface displacement sig¬ 
nals of interest appear as either control or error signals which may be easily picked 
off of the feedback control circuit. The mirror drive control signal follows and may 
therefore be used as a measure of low frequency surface displacements (Leiner, 
1978). The feedback error signal, on the other hand, tracks pathlength variations 
resulting from acoustic waves, for example, which occur at frequencies beyond the 
response limit of the feedback system (Palmer, 1977; Palmer, 1979; Smeets, 1977; 
Emmony, 1982). Techniques other than feedback stabilization have been used to 
ensure peak sensitivity of the interferometer. In one system developed at RCA 
labs (Mezrich, 1974) the reference mirror is moved at a low frequency and rela¬ 
tively large amplitude to assure that at some point in the measurement period the 
interferometer is in its optimum configuration. This technique works espesially 
well for measuring continuous, periodic signals. In such cases the largest signal 
recorded during the test interval corresponds to the point of maximum sensitivity 
and is retained for analysis. Another technique introduces a 90° phase shift 
between orthogonal polarizations of the reference beam. When the object beam is 
interfered with each of the reference components, the two resulting interference 
intensity signals may be squared and then added thus linearizing the interferometer 
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at its peak sensitivity (Vilkomerson, 1976; Peck, 19S3). Finally, optical hetero¬ 
dyne techniques may be used to minimize sensitivity variations caused by path- 
length instabilities. This technique will be reviewed in some detail in following 
sections. 


It is usually assumed that in the case of the Michdson interferometer the shot 
noise of the photodetector sets the fundamental limit to the interferometer’s sensi¬ 
tivity (Kwaaitaal, 1980). This may not always be the case, however. Other contri¬ 
butions to the total noise of the interferometer system include laser noise, en¬ 
vironmental (mechanical) noise, and electronic noise. Kwaaitaal’s work indicates 
that laser noise can be made insignificant if a single mode laser is used and the fre¬ 
quency content of the detected signal is above 10kHz (Kwaaitaal, 1980). Electron¬ 
ic and environmental noise, however, are not so easily neglected. 


Consider first the noise associated with a photodetector/amplifier combination 
required to obtain a voltage signal analog to the light intensity incident on the 
detector. Two general classes of photodetectors will be considered - photoemissive 
and solid state. Specific representatives of these two classes are the photomulti¬ 
plier tube and the junction photodiode, respectively. Tn both cases light incident 
on the active regions of these devices (i.e. photocathode or depletion region) li¬ 
berates electrons or generates electron-hole pairs which contribute to current flow 
through the device. This signal current may be described as follows: 


' h 


(5) 
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• 

where 


£-• 

if — quantum efficiency in 

■- , • -i 
-;o 

m 

carhers/phoion 


q — electron charge 

• 


P — incident light intensity 





■ 

v — light frequency. 

Owing to the quantum nature of light there is an associated shot noise with a spec- 
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iral density function 


■ 

S^y-lNq 1 , (6) 






where N represents the average carrier emission rate (Yariv, 1976). Since the 

average current/, is qN, the shot noise power for a given frequency bandwidth, H, 


B 

ft 


can be written: 



ij-lqlB. (7) 


Lf 


B,.,] 


Thus if shot noise were the only source of noise in the detector/amplifier system. 


• 

then the signal to noise ratio could be computed from Equations 5 and 7 as fol¬ 

• 


lows: 


• 

' 2*(t wP/hv)B 

S/N ~£a' (,) 


• 


•_ 







The quantum efficiency of common photocathode material* is typically about 10% 
versus about 30% for a silicon photodiode (Anderson and McMuruy, 1966). Thus 
if one considers that the performance of a detector/amplifier system is shot noise 
limited, then better signal to noise ratio will be obtained using solid state photo¬ 
diode detectors. 

Tn a shot noise limited system the minimum detectable surface displacement 
may be computed for the case of the Michelson interferometer. Equation 4 pro¬ 
vides a relation for the power incident on a photodetector for a mirror displace¬ 
ment, D. 


P-Pdtl+Mcos ~-D), 

where M (modulation factor) — 2 
PsubO - A^+A^ 2 . 

In the case where an acoustic wave is launched on the mirror surface, the displace¬ 
ment, D is the sum of a static component and a dynamic component. Using 
Kwaaitaal’s notation; 


D-X+x 


where X — mirror position 


x —isin 2irft. 
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Thus the power incident on the photodetector is 
P—/ 3 0 (l-hMcaj[ ^Uf+*)]) 

-Pdtt+Mco3{~-X) cos^x>-Msin{~X) sin(~x)) 


For x<<X, 


P-PJLl+Mcasi^-X )-M~-ain {~-X )). 


(9) 


Note that only the third terra contains the signal displacement factor, x. Conse¬ 
quently the average photocurrent contributing to detector shot noise derives from 
the first two terms of Equation 9: 


ij$—2qIB 

-Iqi^qPjh „)( \+Mcos{~~X))B. 


The signal current is 


i,—(^ Wll^-X ). 

hv X X 


The signal to noise ratio is now computed as 


SM 


I# X 


x) 1 - 




l-hWcoj^Ar 


( 10 ) 
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For other than ideal systems, M< 1 causing the maximum signal to noise to occur 
for jr=3X/16. This result is shown graphically in Figure 2a. Note that, somewhat 
surprisingly, the best signal to noise ratio is not obtained when x— A/8, i-e. where 
the slope of the signal curve is greatest (Figure 2b). Instead the best SNR is ob¬ 
tained for A"=3X/16. One can reconcile this result by observing that the 3A/16 
point represents a compromise between maximum signal gain and minimum shot 
noise. The background light level at the point of maximum signal slope (X—A/8) 
is greater than at the computed SNR maximum. Therefore the average photodetec¬ 
tor current and consequently the shot noise will be greater at this point. On the 
other hand, when X—A/4, the shot noise is at a minimum. Unfortunately the sig¬ 
nal gain is at a minimum as well. Therefore the optimum value for X must lie 
somewhere between A/8 and A/4. Using the argon ion laser line at 514.5nm, and 
a 0.5mm diameter silicon detector, and 0.05mW/cm 2 output field intensity typical 
of a holographic interferometry system, the minimum dynamic displacement sensi¬ 
tivity for a SNR—1 can be computed using Equation 10. For X—3A/16, M—0.9, 
and a 1Hz bandwidth. 


x—4.71xl0“ 8 m. 

Note that this is about 500 a and is, therefore, far less sensitive than desirable. 
Fortunately, sensitivity may be improved by increasing laser illumination, holo¬ 
gram efficiency, detector area, and by using heterodyne detection all of which tend 
to increase the light intensity incident on the photodetector. Bear in mind that this 
rather disappointingly poor sensitivity is largely the result of the low light intensity 









figure which was selected to represent actual intensities encountered in holographic 
interferometry. 


Heterodyne Interferometry 


Heterodyne interferometry provides another alternative to path stabilization for 
immunity to environmental noise sources. In addition, the application of hetero¬ 
dyne techniques to holographic interferometry allows for interpolation between 
fringes to obtain sensitivity to static displacements comparable to real time dynam¬ 
ic detection sensitivities obtained by classical interferometers. The basic difference 
between heterodyne and homodyne interferometry is that the two wave fronts to 
be interfered in the output plane of the interferometer are shifted in frequency re¬ 
lative to each other for heterodyne and not for homodyne. The light intensity in¬ 
cident on the photodetector may be computed in a manner similar to that used in 
the development of Equation 3. 


Ut-A* 1 '’***' 


Then the intensity, T, at the detector is: 


HUi+U jMU'i+U'jXU ,-W 2 ) 







•A f+A^+lA \A^03 [(o>,— 


( 11 ) 


The photodetector signal, then, is an offset sinusoid whose intensity varies at the 
difference frequency, /—27r(a»,—*» 2 ). When used in a Michelson interferometer 
configuration, a dynamic displacement of the object mirror results in a phase shift 
in the detected signal. Therefore, displacements are detected using electronic 
phase measurement techniques which can be made insensitive to signal amplitude 
variations. 

Several techniques have been used to detect and process the phase shifted pho¬ 
todetector signal. Frequency demodulation has been used to measure the 
Doppler-shifted output signal (Whitman, Laub, and Hates, 1968). The doppler 
shifted frequency is proportional to the velocity of the surface motion and there¬ 
fore electronic integration must be performed in order to obtain displacement in¬ 
formation. The difficulties in accurately performing the required integration partic¬ 
ularly at frequencies very near the difference frequency make Doppler frequency 
detection a less attractive option for heterodyne signal processing (Palmer and 
Green, 1979). FM detection using a phase locked loop has been used to success¬ 
fully detect high frequency vibration and pulsed displacements (Rudd, 1980). In 
such a system an electronically phase locked loop is locked to the heterodyne sig¬ 
nal from the photodetector. Phase variations occurring at a rate below the loop 
response rate (typically 100kHz) are tracked by the loop. Thus the output of the 
phase detector remains constant for low frequency phase changes. Since the fre¬ 
quency content of most environmental noise sources is well below 100kHz, the 
loop provides a level of noise immunity si milar to that obtained with path length 
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stablilized interferometers. For vibrational disturbances above 100kHz, the error 
signal present at the output of the phase detector provides a voltage signal propor¬ 
tional to the velocity of the vibrational disturbance. Another scheme also sensitive 
to the velocity of the test object's surface displacement uses coherent detection of 
the phase shifted signal (DeLaRue, 1972). This technique works particularly well 
when the frequency of the acoustic exciting source is well known. In this case the 
signal to be detected is shifted to a sideband of the detector carrier frequency. By 
using the earner signal as a local oscillator, immunity to environmentally induced 
path length instability is obtained. This optically derived local oscillator is then 
mixed with one of the signal sidebands to produce the desired output signal. In 
another coherent detection scheme for phase demodulation, the local oscillator is 
simply the input signal delayed by a known amount, typically = 500ns (Jacobson, 
1982). By delaying the local oscillator by some integral number of wavelengths 
plus A/4, it can be made to be in quadrature with the heterodyne signal. Coherent 
detection using these two signals and subsequent filtering can be used to to provide 
a direct measure of phase. Again low frequency noise due to pathlength instability 
may be filtered out. 


When changes in the detected phase difference occur slowly, direct phase 
measurement of the detected carrier frequency may be used. To do so requires a 
local oscillator source coherent with the detector signal. If the relative phase 
difference between the two is measured, then any change in the optical path 
difference between the reference and object beams will cause a corresponding 
change in the measured phase difference. Electronic phase difference measure¬ 
ments are quite easily and most commonly made using electronic instrumentation 
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which produces an analog output voltage proportional to the time between 
corresponding zero crossings of the two input signals expressed as a fraction of the 
signal period. This type of phase meter is capable of providing accuracies better 
than 0.1° (Mastner, 1980) thus giving a potential optical path difference sensitivity 
of (A/2X0.1/360) — X x 1.4 x 10" 4 or about 0.07nm. For the work reported here, 
the relative phase differences between the signal and reference oscillator were 
derived from the outputs of a quadrature detection lock-in amplifier. Using a 
lock-in amplifier also reduces the effects of variations in input signal amplitude en¬ 
countered with other phase measurement schemes (Ineichen, 1977). 

Coherent quadrature detection for relative phase measurements is shown 
schematically in Figure 3. The input signal is mixed both with the local oscillator 
and with a 90° phase shifted version of the local oscillator. 

Signal, S — AcosCwt 4- 0) + N 

Local oscillator, R ( — He ascot 
Shifted L.O., R 2 — Hsirvut. 

Output from first mixer is 

O \—SR ^ABcosiotcos{o>t-H>}i-NBcostot 

—( )cos (2wf -h£ >-K )cos dH-NBcos ctf. 


Output from the second mixer is 











Acos(wt+<£)+ N 
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Oi-SRfABsiniut cos(**+6h-NBsin*it 
-i ^ hin ( 2atf +* W ^ )rm <H-NBsn <w. 


Note that by low pass filtering the outputs O ( and O 2 only the second term in 
each expression is retained (assuming that A { and.4 are slowly varying func¬ 
tions of time). Note also that the additive noise, N, originally present in the input 
signal is eliminated by this filtering process as well. In fact it is by this detection 
and filtering process that lock-in amplifiers may detect and amplify signals in noise 
backgrounds as much as 1000 times greater than the signal amplitude! The signals 
remaining after the filtering process then are as follows: 

OF x -~cosA 

OFi-^arxf, 


By taking the arctan of {OFjOFi), the relative phase difference, <t>, is computed. 
Some lock-in amplifiers are available which compute arctan (OF JOF ,) by an ana¬ 
log process in real time. Again owing to the coherent detection and filtering pro¬ 
cess, the ratio ( OFjOF \) and <£, therefore, is independent of the signal amplitude, 
A. 


Judicious choice of the local oscillator source can reduce the effects of environ¬ 
mental noise on the accuracy of the relative phase measurement. The most com¬ 
monly used technique for providing optical frequency shifting in heterodyne inter¬ 
ferometry is the use of acousto-optic modulators, or Bragg cells. Other techniques 
which have been used include rotating diffraction gratings, routing wave plates. 
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and electro-optical "serrodyning" (Hard, 1978; Aleksoff, 1971). Using a Bragg cell, 
the first order diffracted beams exiting the cell are shifted into upper and lower 
sidebands around the optical frequency a>. Light in either of the diffracted orders 
varies in frequency from the undiffracted beam by ±/ D , where f D is the RF driv¬ 
ing frequency of the cell. The frequency shift is the result of a phenomenon 
known as the Debye-Sears effect (Adler, 1967). Typically the RF drive frequency 
for commercial Bragg cells is 40 or 80MHz. Therefore if an 80MHz Bragg cell is 
used to perform the frequency shifiting in a heterodyne interferometer, the 
difference frequency signal from the photodetector will be equal to the 80MHz 
drive frequency. While a high frequency carrier may be desirable in real-time 
detection of acoustic displacements, it is not necessary for the detection of slowly 
varying changes in optical path length. Indeed the design and construction of 
amplifiers and associated processing electronics becomes somewhat more compli¬ 
cated in the higher frequency regimes. 

A more manageable carrier frequency is obtained by passing the laser beam 
through two Bragg cell modulators driven at slightly different frequencies as shown 
in Figure 4. Light diffracted into the +1 order of the first cell is frequency shifted 
upward by 80MHz. This first order diffracted beam passes subsequently through 
the second Bragg cell. The -1 order beam exiting the second cell is now shifted 
down in frequency by 80.1MHz. Therefore the net frequency shift caused by pass¬ 
ing the beam through both cells is only 0.1MHz (100kHz). The resulting 100kHz 
difference frequency detected at the output of the heterodyne interferometer is far 
easier to amplify and process than would be the 80MHz. The diffraction efficiency 
of commercially available Bragg cell modulators using lithium niobate is better 
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Figure 4: Obtaining a 100kHz frequency 
shift using two Bragg cells 
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th&n 80%. Consequently, the overall efficiency of this frequency shifting process is 
on the order of 65%. Unfortunately it is impractical to use just a single cell at 
100kHz driving frequency because the diffraction efficiency would be quite poor, 
and more importantly, an extremely small diffraction angle of the 100kHz shifted 
components makes it difficult to spatially resolve the diffracted and undiffracted 
orders. There is, however, one very convenient feature of practical importance as¬ 
sociated with using two Bragg cells for frequency shifting. That is the fact that both 
cells may be driven at the same shift frequency resulting in no net frequency shift 
nor causing any appreciable change in the optical path or angle through ’he cells. 
Therefore the heterodyne interferometer may be used in a homodyne mode 
without the need for reallignment of the beam in the shifted leg. This feature is 
particularly useful in initial aliignment of the heterodyne interferometer. 

Most of the discussion to this point has been with regard to the use of inter¬ 
ferometers to detect and measure dynamic displacements at a single point on the 
surface of an object mirror. In fact, the research reported in this paper also ad¬ 
dresses the problem of detecting and mapping dynamic surface displacements asso¬ 
ciated with traveling surface acoustic waves. However, unlike the imerferometric 
techniques which have been discussed to this point, heterodyne holographic inter¬ 
ferometry is not limited to real time surface displacement measurements at a single 
point on the test object. Instead, the surface displacements over the entire field of 
view of the test object are recorded at a single moment in time using pulsed holog¬ 
raphy. Subsequently, heterodyne interferometric techniques are applied to the 
reconstructed holographic image to provide details of the surface deformation over 
the entire field of view. This later step is quite similar to broad beam heterodyne 





interferometry which has been used to measure surface flatness of optical com¬ 
ponents, for example. 

In broad beam interferometry, the illuminating beams are expanded and col¬ 
limated so that large surface areas may be Lnterferometricaily compared. If the ob¬ 
ject minor is tilted or deformed relative to the reference mirror, a fringe pattern is 
generated in the output plane as shown in Figure 5. The fringes map out contours 
of constant displacement with a contour interval of A/2. When a relative frequen¬ 
cy shift is introduced between the reference and object beams, a detector placed in 
the output field generates a sinusoidal signal at the shift frequency. Direct meas¬ 
urement of the phase difference between this heterodyne signal arxi a local refer¬ 
ence oscillator may then be used to interpolate between fringes providing a direct 
measure of out of plane displacement as a function of position on the object mirror 
surface. By scanning the detector over the output field, a map of surface topogra¬ 
phy may be generated. Such real-time broad beam, heterodyne interferometry has 
been used to measure flatness and curvaiure in mirrors, silicon wafers, and optical 
components (Massie, 1979). 

Broad beam interferometry provides yet another possibility for reducing noise 
resulting from path instabilities in interferometry. A reference detector may be 
placed at any fixed point in the output field. The sinusoidal signal from this refer¬ 
ence detector may be used as a local oscillator signal against which the phase of the 
signal from the scanning detector may be measured. Low, frequency vibrational 
disturbances which may affect the optical pathlengths of either the reference or ob¬ 
ject leg of the interferometer will cause phase shifts in both the reference and 
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scanning detector signals. Since both detectors are placed in the same output 
plane, the phase shifts resulting from vibration will be equal in both signals. 
Therefore the easured phase difference between the two signals will be 
unaffected by the vibrational disturbance. 


Computing the noise limited sensitivity of a broad beam heterodyne inter¬ 
ferometer requires an understanding of the effect of signal to noise ratio (SNR) on 
the performance of the system’s phase meter. Assuming that the detector 
bandwidths are less than the heterodyne carrier frequency, then a phase error, A£, 
using a zero-crossing phase meter is 


where r— integration time 
f — heterodyne frequency 


(Ineichen, 1977). 


If phase fluctuations from vibration and other environmental effects occur at fre¬ 
quencies outside of the system bandwidth, then the SNR can be assumed to 


depend only upon shot noise. 


To compute the shot noise limited SNR, recall from Equation 11 that the light 
intensity incident on the photodetector is 


P-P 0 (l+Mcwt(» r -ai 2 )r4*]) 
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The phase error resulting from a shot noise limited system which incorporates a 



zero-crossing detector phase meter is 
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Using the same values for tj, M, P 0 , and B as used earlier for the example calcula¬ 
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tion for the path-stabilized interferometer (i.e. 17 — 0.3, M ” 0.9, Po*“10 ~ 7 watts, 



and B — 1Hz) and using an Ar* laser (X — 0.5145micron) and a zero-crossing 
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phase meter where /r—l, the phase error, A<£, for a heterodyne interferometer is 
8X10 -6 radians corresponding to a displacement sensitivity of about 0.003 a. Note 
that this is some 5 orders of magnitude better than was obtained for the path sta¬ 
bilized interferometer. 


Holographic Interferometry 


To this point, all discussion of interferometry, path stabilization, and detection 
has been illustrated by consideration of the Michelson interferometer. Holographic 
interferometry incorporates many of the same operating principles but reduces 
some restrictions imposed by classical interferometry, particularly regarding the 
types of specimens which may be examined. For any form of interferometry to be 
of practical use, the interference fringe pattern in the output plane must be of rela¬ 
tively low order so that individual fringes may be resolved. Obviously if the fringe 
rate is too high, there will be no fringe pattern observed visually or detected with a 
photodetector of finite area with dimensions larger than one fringe period. This 
necessity for a low fringe rate imposes the restriction that the two wavefronts to be 
interfered must be very nearly identical. Tn the case of a classical interferometer, 
nominally identical wavefronts are reflected from two flat mirrors with one mirror 
serving as a reference. When the object mirror is moved, tilted, or deformed, 
differences between the reference and object wavefronts occur which manifest 
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thernselves in the observed interference pattern. If the object mirror is replaced by 
a diffusely reflecting surface, a very complicated and high fringe rate pattern ap¬ 
pears in the output plane. The source of this complex pattern is clear if one con¬ 
siders that a diffusely reflecting surface is made up of a mosaic of small mirror-flat 
portions at random angles each reflecting a portion of the illuminating beam back 
to the output plane. Therefore the fringe pattern in the output plane is the result 
of interference of the light reflected from each of the individual reflectors with the 
reference beam as well as with each other. (It is the interference pattern from the 
interaction of the light from the individual.surface reflectors with each other which 
is called speckle. The effects of speckle on the accuracy of holographic inter¬ 
ferometry measurements will be discussed in some detail shortly.) By taking ad¬ 
vantage of the "mosaic" nature of diffuse surfaces, one may still use classical inter¬ 
ferometry to transduce surface displacements. This is often dons by focusing the 
object illuminating beam down to a small point such that at the point of illumina¬ 
tion the local variation of the object surface is negligible. The light is reflected 
directly back from the test point and recollimated by the focusing objective lens. 
Subsequent interference of this object beam with the reference takes place as 
though the beam were reflected from a flat mirror. 

There are a number of applications, however, where it is desirable to inter- 
feroraetrically examine an entire surface. As discussed previously, broad beam in¬ 
terferometry may be used for objects with specularly reflecting surfaces which are 
either flat or of a known simple geometry. In such a system, a "master" object 
against which a test object is to be compared is used in place of a reference mirror. 
Usually wavefront correcting lenses are placed in each leg of the interferometer so 
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that the extent of the beams reflecting from each object can be controlled to stay 
within the physical dimensions of the components comprising the interferometer. 
The use of conventional spherical lenses restricts the nominal shapes of the test 
object surfaces to be either flat, convex, or concave. 

Interferometric analysis of more complex, diffusely-reflecting surfaces requires 
that a reference wavefront be generated which matches that of the test object sur¬ 
face One might consider the possibility of manufacturing a roaster object to re¬ 
place the reference reflector in a Michelson interferometer as described above It 
is important to remember, however, that in order for such an approach to work, 
every dimension of the master object including each minute reflector in the "mosa¬ 
ic* comprising the object surface must match that of the test object to be exam¬ 
ined. It is certainly impractical, if not impossible, to produce reference surfaces 
with the required precision. 

Holography provides one solution to the problem of generating a master wave- 
front for interferometric comparison of objects with arbitrary surface shape and 
texture. While some work has been done adapting holography for interferometric 
comparison of different, nominally identical objects (Neumann, 1981), by far the 
most common applications of holography in interferometry is in real-time or 
double-exposure interferometry where the object surface is compared at two 
different times. By interfering the wavefronts which reflect from the object both in 
a stressed and an unstressed state, for example, one may observe deformations 
which occur in the object as a result of applied stress. If the image of an object as 
recorded on a hologram is interfered with the image of the object at a slightly 






different longitudinal magnification, the object appears superimposed with interfer¬ 
ence fringes forming contours of constant surface elevation. 


The operating principles of holography may be explained using either wave- 
front interference theory from geometrical optics or alternatively using scalar 
diffraction theory. It is the latter approach which will be pursued here. In scalar 
diffraction theory, often called Fourier Optics (Goodman, 1968), an optical wave- 
front may be described as a complex quantity whose magnitude and phase varies as 
a function of x and y in a plane normal to the optical axix, z. A phasor notation of 
sorts may be used to explicitly specify the propagation direction of a wavefront as 
the angle between the propagation direction and the optical (z) axis as follows 
(Collier, 1971): 

where 4*-/(0)a0. 

Looking at the optical set up shown in Figure 6, note that two wavefronts intersect 
at the film plate. One wavefront results from the reflection of the illuminating 
beam by the object surface and will be referred to as the object wavefront, 0 . The 
second beam is unmodified by reflection and is called the reference beam, R. 
Defining the optical axis to be perpendicular to the film plate: 


R~Re J * 

6 -Oe°-0. 
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Thus the optical field falling on the film plate, A —{&*+0. Assuming that the 
transmittance, T, of the exposed and subsequently developed film is proportional 
to the field intensity, I, then: 

-(/S+OX/T-K?*) 

H£e J *-K? XR'e~ J *+0 ') 

-*R\ l MO\ l -OR' e~ J *+Q' R* J *. (13) 

Upon reillumination by the reference beam as shown in Figure 7, the hologram 
diffracts the incident light into three beams corresponding to the three bracketed 
terms below. 

Tx ; «e>M(| ; R| z +|0| 2 )^*}HO| : R| 2 H-[0*l^|V*] (14) 


Analyzing Equation 14, the first term describes the intensity variations made to the 
reconstructing reference beam as it passes through the hologram plate and contin¬ 
ues to propagate at an angle, 9, to the optical axis. The third term contains the 
conjugate image of the object and leaves the film plate at an angle 29 relative to 
the optical axis. It is the second term of Equation 14 that is of importance in holo¬ 
graphic interferometry. Note that this term describes a wavefront which pro¬ 
pagates along the optical axis (<C*- 0 — 0 ) and is the product of the original object 
wavefront multiplied by the squared magnitude of the reference wavefront. In the 
case being considered here, the reference wave is a uniform, collimated beam such 
that its squared magnitude is simply a scalar constant independent of x and y. Thus 
a reconstructed replica of the original object wavefront propagates along the optical 
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axu such that an observer portioned as shewn in Figure 7 would see a virtual im¬ 
age of the object appearing in its onginal location and in 3-D. 

The hokigraphscaiiy recorded image of the object may new be used as a 'mas¬ 
ter' in an interferometer as shown in Figure 8. The recording and playback 
geometries shewn in Figures 6, 7, and 8 were chosen to illustrate the similarity of 
holographic interferometry to the classical Mi c he Ison interferometer. Tn practice, 
however, off-axis illumination of the test object is used and the hologram serves as 
both a beam splitter and reference image source. Off-axis illumination does, how¬ 
ever, complicate fringe analysis somewhat. Remember that interferometric fringe 
formation is the result of optical pathlength differences which occur between the 
reference and object wavefronts. For coaxial illumination, i.e. when the illumina¬ 
tion and viewing axes are the same, one fringe corresponds to an object displace¬ 
ment of A/2. For off-axis illumination, however, the sensitivity to object displace¬ 
ment is somewhat reduced. For analytical purposes a sensitivity vector, K, is 
defined in the direction of the bisector of the illuminating and viewing axes (see 
Figure 9) (Vest, 1979). The sensitivity vector has a magnitude; 

[K|-2Aroj0 

where k—-— 

0 — half angle between the 
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Figure 9: Sensitivity vector, K, for off-axis 
holographic interferometry. 
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For an object displacement D, the observed intensity in the output plane of the 
holographic interferometer varies in a manner similar to that described in Equation 
4 as follows: 

I-2A\l+co3 K D ) (15) 

—lA^l+cosOkDcas flcosd)], 

where £ is the angle between the sensitivity vector and. the displacement direction, 
D. Another advantage of holography is now apparent. Using Equation IS the 
component of D in the K direction may be determined for any point on the object 
surface. By selecting three different viewing angles, the projection of D on three 
different K vectors may be computed and subsequently processed to completely 
specify D in magnitude and direction. 

Double exposure holography eliminates the need to precisely replace and align 
the hologram in the optical set-up after its exposure and development. In this case, 
both the reference and object images are reconstructed by the hologram. The opti¬ 
cal arrangement for recording a double exposure hologram can be identical to that 
used in real lime holographic interferometry. Two exposures of the object image 
are recorded with a change either in magnification or applied stress being made 
between exposures. The alignment of the reference and object images is preserved 
in the hologram so that upon reconstruction with a single reference beam, both 
images are superimposed and therefore interfere. One variation of this double- 
exposure geometry is important in heterodyne applications. This variation is 
shown in Figure 10 where two spatially separated reference beams are used to 
record and reconstruct the holographic images. One reference beam and then the 
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Figure 10: Dual reference beam holography 
for heterodyne interferometry. 









other is used in making the two exposures respectively. On playback, both beams 
are used to simultaneously reconstruct the reference and object images. If one of 
the reconstrudng reference beams is shifted in frequency relative to the other, 
then heterodyne analysis techniques similar to those used in broad beam hetero¬ 
dyne interferometry may be used. 


Heterodyne Holographic Interferometry 


Heterodyne analysis of dual reference beam, double exposure holograms 
recorded as described above, permits interpolation between the normally observed 
interferometric fringes. By interpolating between fringes both the sensitivity and 
dynamic range of holographic interferometry may be greatly improved- When first 
introduced in 1973, interpolation to better than 1/100 of a fringe was demonstrat¬ 
ed (Dandliker, 1973). In 1980, Dandliker published a rather complete review pa¬ 
per on heterodyne holographic interferometry which predicts practical sensitivity 
limits of better than 1/1000 of a fringe (Dandliker, 1980). Indeed Dandliker and 
his associates dominate the literature from 1973 to 1980 on this topic. Consequent¬ 
ly much of the theoretical development to be presented here is discussed in greater 
detail in his 1980 publication. 


In many respects, the principles of heterodyne holographic interferometry 
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paiaUel those of broad beam heterodyne interferometry discussed earlier. In the 
holographic case, however, the two wavefronts to be interfered are reconstructed 
from the hologram plate rather than being reflected from reference and object mir¬ 
ror surfaces. By constructing the hologram with two reference beams as shown in 
Figure 10, the two reconstructed image wavefronts interfere only when both refer¬ 
ence beams illuminate the developed hologram simultaneously. When one refer¬ 
ence beam is shifted in frequency relative to the other, the two reconstructed im¬ 
age wavefronts are likewise shifted in frequency causing a corresponding variation 
in the image intensity. Mathematically, the transmittance of the dual reference 
beam, double exposure hologram can be expressed as the superposition of the 
transmittances which would result from each exposure independently as shown 
previously in Equation 13. 

Ta (16) 

where the subscripts, 1 and 2, refer to the first and second exposures, respectively. 
Again, only the last two terms will give rise to wavefronts which will propagate 
along the optical axis to the output plane. 

To perform heterodyne analysis for fringe interpolation, the two reference 
beams used to reconstruct the object wavefront, O, and 0 2 , must be at slightly 
different frequencies. Let 


R,-R 
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When together they illuminate the hologram whose transmittance is as shewn in 
Equation 16, the light intensity falling on the ouiput plane of the interferometer is 
obtained as follows: 

H(/?i+K 2 )x7’| 1 

+0J/?JV 3< 

+o I l* 1 |V ,, | J . 

If R , and R 2 are plane waves of unit amplitude then 

HO J /"i»**i**>*<7 

Note that only the last two terms of this expression correspond: to wavefronts 
which propagate exactly down the optical axis toward the output plane. The first 
two terms form images which lie to either side of the optical axis propagating at an 
angle equal to that between the reference beams. Therefore considering only the 
last two terms above and noting that 

Of-o*'****'. 


Then 
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“0?+02+20i0 2 c<M[(<ui-^a» 2 )<-l-^ 1 —^j]. (17) 

This is exactly the same form as obtained in Equation 11. Therefore sub-fringe in¬ 
terpolation between interferometric fringes may be accomplished at any point in 
the image by scanning a detector over the image plane and applying any of the 
phase detection schemes discussed previously. 

Unlike the case of simple heterodyne interferometry, the resolution of hetero¬ 
dyne holographic interferometry is limited both by shot noise and speckle noise. If 
a zero-crossing phase meter is used to measure the phase difference between sig¬ 
nals from two detectors located in the output plane, then the shot noise limited 
phase sensitivity is the same as that of broad beam heterodyne interferometry as 
expressed in Equation 12. To determine the effect of speckle on the accuracy of 
phase measurements it is helpful to understand some of the statistical properties of 
speckle The discussion which follows is not intended to provide a complete 
description of speckle phenomena but rather will attempt to develop those proper¬ 
ties of speckle which directly affect the accuracy of heterodyne holographic inter¬ 
ferometry. 

Speckle is the name given to the phenomenon which gives a granular appear¬ 
ance to diffusely reflecting (or transmitting) objects illuminated by a coherent 
source. The origin of speckle is quite easily explained when one considers that a 
diffusely reflecting (rough) surface can be considered to be composed of a mosaic 
of mirror-flat sections whose dimensions are on the order of the wavelength of the 





illuminating Light and whose faces are oriented randomly about the average surface 
normal. At any point in front of such an object, light from the illuminating source 
is scattered by a number of surface reflectors creating a complex interference pat- 
terra In general, this speckle pattern is observed using an imaging system such as 
the eye or a camera so that there is a finite aperture through which the light must 
pass before it is focused on the image plane. This finite aperture limits the resolu¬ 
tion of the imaging system thus establishing the speckle size at that limit. Conse¬ 
quently for a square or a circular aperture, the speckle size is Xz/d or 1.22Xz/D, 
respectively, where z is the distance from the aperture to the image plane and D is 
the aperture width (Jenkins and White, 1950). The results of this rather intuitive 
argument for determining speckle size agree with those obtained by more rigorous 
analysis (Dandliker, 1980; Vest, 1979; Goodman, 1975). Speckle size, as will be 
shown shortly, is an important factor in determining fringe contrast of holographic 
interferograms. 

As discussed earlier, the only practical way to interfere wavefronts reflecting 
from diffuse object surfaces is through holographic interferometry. In the 
mathematical description of the two reflected fields, a roughness term, p, must be 
added to account for the microscopic variations in surface topography which give 
rise to speckle. Therefore the two wave fields are expressed as follows: 

^2"«2^o)Pl(io). 


where and a 2 represent the macroscopic, or smooth, variations in the reflected 
wavefront and xq—( xo^o) positions in the object plane. On the output plane of a 
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iinear space-invarient imaging system, the complex amplitude of each wavefront 
may be computed by convolution with the impulse response of the imaging sys¬ 
tem, h(x ): 


£ t (x,>“fa ](x 0 )pi(x 0 )/i (x, -Mx 0 )dx 0 (18) 

B 1 (x, >-J*i a 2 (x 0 )pj(xo)/t (x, -A6c 0 )dx 0l 


where i>Kx,»y,) positions in the image, or output, plane of the imaging system 
and M is the lateral magnification of the system. The intensity in the image plane 
is the coherent superposition of the wavefronts 5, and B 2 such that 

HB x +B ,1 x B\+B jBj. 


As always, the last two terms represent the interference of the two wavefronts. To 
separate the microscopic (speckle) interference from the macroscopic interference 
it is useful to examine the average interference field over a large number of 
diffusers, p. The emsemble average of the interference terms may be expressed as 

<5 \Bi~\~BjB i >. 

However, since <B X B\> — <B^\ > , the overall average will be twice the aver¬ 
age of either term. Therefore consider 


<B >—< ( ft a i(xo)pi(£oV* (x )— A/xo)<&o )x( J* o 2^io)pj(xo)** (x, — Mx 0 Vx 0 )>- 



By changing the dummy variable in the second integral to xo', the order of multi¬ 
plication and integration may be reversed such that 




Finally, p\ and p 2 are considered n be ergodic random processes which is to say 
that the statistical properties of th roughness over any sample area of the object 
surface are equal to those obtained for any other sample area. Consequently it is 
permissible to reverse the order w :h which the integration and averaging is per¬ 
formed (Lathi, 1968) so that 

<B\B\>—ffa |(xo)aj(x 0 ')<Pi( )pj(x 0 ')>/i(x, —A/x 0 )/»* (x J —Mx 0 ')dx 0 'dxoi 19) 

Through Equation 19 it is possii i to demonstrate a great deal about the effect 
of surface roughness and speckle n the interference of two diffusely reflected 
wavefronts. First, to convince oneself of the physical significance of Equation 19, 
consider the case where B \=-B 2 , then 

<B\B\ > mm <B l B’ l >-<lBj 2 >, 

which is the average intensity of B t in the output plane of the imaging system. 
Since Bf-B 2 then Pi—p 2 . Thus <Pi(xo)p*(xo)> is the autocorrelation of the sur¬ 
face roughness R p (x o~xo’). Recall that the spatial frequency spectrum of the sur¬ 
face roughness is assumed to extend well beyond the bandwidth of the imaging 
system. In other words the fine detail of the surface cannot be resolved by the im- 








aging system. Therefore the autocorrelation peak of the surface roughness will be 
much narrower than the resolution of the imaging system and can be approximated 
by a delta function: 


<Pi(*o)piQto ) >= ^*cr2o)- 

Owing to the shifting property of the delta function, Equation 19 is considerably 
simplified, 

<|5| 2 >-J*J*a,Gf 0 )ai (x 0 )Mx;—Mxo) 

x h'ix,-Mx 0 )b{xQ—XQ)dxtfix 0 

-/1 a il HxoM -Mxo'X&o'- 

The above expression is now recognized as the convolution of the input intensity, 
|a| 2 , with the incoherent transfer function. 

Having established the physical significance of Equation 19, return now to the 
interference case where B^B V First observe that, for interference to exist at all, 
the roughness terms for the two wavefronts must be correlated, i.e. 
<PiCxo)P 2 (io ')>^0 • practice, then, one can interfere two speckle fields only 
if both wavefronts originate by reflection from the same object surface in two 
different states. This explains the value of holography in performing inter¬ 
ferometric analysis of diffuse object surfaces. Assume that following the recording 
and subsequent development of a hologram, the test object is translated by an 
amount u as shown in the real time holographic interferometer of Figure 11. The 
interference terra (Equation 19) will be affected in two ways. First, Bi will be 
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Figure 11: Real time holographic interferometer. 
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shifted in phase by an amount equal to the projection of u on the sensitivity vec¬ 
tor, K (see Figure 9). Secondly, there will be a lateral shift in the speckle field of 
B 2 in the image plane by an amount equal to the product of the lateral 
magnification, M, and the projection of u onx<> Thus 

5 2 (£ >~B ,(x;+M(u •x 0 ))e' ( * !4) 

where 

u, —M (u x 0 )<fr-K u. 

Substituting into Equation 18: 

Bf-e j4 fa iC*o)Pi(*oM (*, +«, -iVfxoVx 0 . 

So that Equation 19 becomes 

<B X B\ i(x 0 )al (xo) <Pi(^o)p*(£o)>^ -Mx 0 )h'(xj+Uj-Mx 0 )dxo'dxo 

—e J *ffa\(x () )a] (x 0 ')/»(x, —A/x 0 )/»* Or,+t^—Afxo)8(x 0 —x 0 )dx 0 'dx 0 
Ha \\Kxo)e J *f hixj—Mxo )h’ix l +u i —Mxo")dxo' 

Ha.RxoV^u,), 

where /?*(«,) is the autocorellation of the coherent impulse response, h. This au¬ 
tocorrelation is also the impulse response function or point spread diameter and 
determines the resolution of the imaging system which is the speckle size. For a 
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circular aperture, the autocorrelation function is an Airy function whose value is a 
maximum for u,-0 and falls to zero when u,—1.22A/D. Therefore the interfer¬ 
ence fringes will be visible only for lateral shifts of the speckle pattern, u ,, which 
are less than the resolution limit of the imaging system. The fringe contrast, y, 
can then be defined; 


y(«,H^(«,)|<1. (20) 

With the development of the concept of fringe contrast it is now possible to 
address the effect of speckle on the accuracy of phase measurements in heterodyne 
holographic interferometry. Since y takes into account effects of the imaging sys¬ 
tem, it is no longer necessary to define the output wavefronts B j and B 2 in terms 
of the input fields and the system transfer function. Instead let 

5|(x,r)—6l(x)V|(x)e J *' ,, (21) 

B 2 (x ,t )—b 2 ( x )v 2 ( x )e J<-jl >) , 


where 6| and b 2 are the average macroscopic image amplitudes and v t and v 2 are 
complex variables related to the object surface roughness and the resulting speckle 
fields. In order for B\ and B 2 to interfere, however it is necessary that v, and v 2 
be somewhat correlated, thus v 2 can be considered to be the sum of two uncorre¬ 
lated fields whose magnitudes are dependent on the contrast, y (Dandliker, 1980): 


v 2 (x >“yv i( 2)« ; *'K l-y 2 )*V(x) 
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<|vj 2 >—<|vj 2 >—1 
<v,v^ 

<VjW* >H) 

<|v,l 4 >-2. 


Substituting into Equation 21: 

B ,(x ,t )—b ,(x )v (x )e Jm * 

BU ,t )-b 2 (x)[yv (x )e"H . 


The interference terra of the intensity: 

B\B 2 o[ y w*e^-Kl-V)Vwk ; ' ,,J * ,,, \ (22) 


where the first terra is the correlated terra and the second is an uncorrelated term 
which will add randomly to the coherent terra introducing a phase error A<£. For a 
detector area large enough only to receive light from two speckles, the vector di¬ 
agram in Figure 12 applies. The net uncorrelated component is the vector sum of 
the random uncorrelated components of each speckle. The correlated components 
add coherently resulting in a net component equal to the sum of the individual 
components from each speckle. One can see then that the phase error, A 4>, will be 
reduced as the number of speckles falling on the detector irxaeases. From the 
figure, the phase error can be expresses as 


sinAtf— 


i u sintfr 
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For small phase errors sinA<£—A0, so the mean square phase error is 


<A* 2 >' 


<|/ 1( | J sinV> 

<I4I J > 


1 <IU 2 > 

2 <|/ C P> 


To compute A^ as a function of the number of speckles, N, on the detector, it is 
necessary to find expressions for <|/ c | 2 > and <|/J 2 >. From Equation 22 

<l4l 2 >“ £ <U/c m > 

N 

“V Z <v « v n v m v m> 

-y 2 OV<|v|>+iVaV-l)<|v| 2 > 2 
-y 2 N(N+ 1); 

<|/„l 2 >- Z </*>•> 

n—m—1 

-V 

—(1—y 1 ) Z <v„w’v m w m > 

—<1—y 2 )Af<|v| 2 ><|w| 2 > 

-K1—y 2 W, 


»nri finally. 


<A* 2 > 


1 <|/J 2 > 

2 <|/ c | 2 > 

1 N(\ —y 2 ) 

2 y 2 JV(N+l) 


(1-y 2 ) 
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The use of two reference beams to record a double exposure hologram for sub¬ 
sequent heterodyne analysis introduces two additional parameters whose effects 
may degrade the accuracy of the phase measurements. These two parameters are 
misalignment of the reconstructing reference beams and wavelength differeras 
between the recording and reconstruction sources. Both are pertinent to this work 
since in the final experiments holograms were recorded and analyzed in two 
separate optical systems. In addition a frequency-doubled NdtYAG laser 
(X—532nm) was used to record the holograms while reconstruction and analysis 
was performed using the 514.5nm line of an argon ion laser. 

Although Dandliker discusses effects of misalignment and wavelength in his 
1980 review paper, the physical significance of these effects is easily lost in the 
mathematics associated with k-space manipulations and scalar diffraction theory. 
For this reason, the following development is offered using geometrical optics 
theory of diffraction gratings. While this approach does not appear in the litera¬ 
ture, it is never-the-Iess appropriate for this background discussion. 

In one sense, a hologram is a spatial-frequency moduated diffraction grating 
whose carrier frequency is determined by the angle between the reference beam 
and the average object beam direction. In the case of a dual reference beam holo- 
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gmra, two independent diffraction gratings are recorded on the hologram - one for 
each exposure. Given that the hologram may be treated as a complex diffraction 
grating, it becomes clear then that any angular misalignment or wavelength change 
will cause a change in the angles of the diffracted beams exiting the hologram. An¬ 
gular changes result in lateral shifts of the images formed in the output plane of 
the interferometer. Such lateral shifts cause a reduction in fringe contrast accord¬ 
ing to Equation 10 which in turn degrades the phase accuracy as shown in Equa¬ 
tion 23. 

To quantitate these effects, consider the hologram geometry shown in Figure 
13 where R\ and R 2 are plane waves. The hologram is drawn as a broken line 
symbolic of its diffraction grating nature. An imaging lens is located as close to 
the hologram plate as possible. Assuming perfect alignment of the reconstructiang 
reference beams, R i and R 2 , and further assuming no change in optical 
wavelength between hologram recording and reconstruction, the reconstructed ob¬ 
ject beams, O x and 0 2 , propagate at an angle 0 r to the hologram (grating) normal 
and are imaged on the output plane where they interfere. The propagation angles 
for R\, O i, R 2 , and 0 2 obey the grating equation for first order diffraction; 

k-a(sinp R +-sinPo ). 


where the 0 angles are measured with respect to the grating normal and a is the 
grating spacing. In the'case diagramed in Figure 13: 


X —a \(sin{0 l —&'Hsin 0,) 


(24) 
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k—a 2 {sin ( H-*fl B ,) 

Consider the case when both reference beams are slightly misaligned such that 
each is displaced angularly from its bisector by an additional amount A8. In this 
case for /?! the change in 9 r may be computed by direct application of Equation 
24: 


sm0,-sin(0,+-A0, W — -sin (P-S)H —sin(0-<6-l-A8))] 
a i a \ 

2cosy (20/+A0, )sin-|(-A«, Hin(M-A« )-sin( 0-8) 

-2cos-^ 20-25-A8)smy(-A8). 


If A8 is small, then 


-A0,c;s0,=-A8coj(0-8) 

A« 

cos 8, 


Similarly for Ri, 


A0; 2 * 


. * cos(P-S) 
-A5---. 

COS 0 ; 


The total angular separation between the two reconstructed object beams is 


Afy^AS/j - A0/ 2 

—A8—— [coa(0—8>+cos(0+8)] 
cos 9, 
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The corresponding translation in the image plane is 

Uj—\8 t z (25) 

—2zA8cos8. 

For a practical system using a circular imaging aperture so that 

2 , 

* mi r D/Xz 1 Xz 

and A—514.5nm, z— 15cm, D—7.5cm, 8—5°, the maximum angular misalignment, 
AS, which will still result in a fringe contrast y-0.7 may be computed since y—0.7 
corresponds to ujZ)/Az-0.5. By substitution, u f is found to be equal to the 
wavelength (5l4.5nm). Solving Equation 25 for AS , the maximum angular 
misalignment is found to be about lCT^degrees! Such fine alignments are achiev¬ 
able with precision optical mounts. However, this result helps explain why fre¬ 
quent realignment is required with changes in the ambient temperature of the la¬ 
boratory. 

Sensitivity to wavelength changes can likewise be analyzed using the diffraction 
grating analogy and Figure 13. This time it will be assumed that ASH) and that a 
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differem wavelength is used for reconstructing the holographic images than was 
used to record the hologram. The reconstruction angle, 9 r , is shifted when the 
hologram is reconstructed using X 2 from what it would have been if reconstructed 
using Xj by an amount which may be computed using Equation 24: 



The same result is obtained for the second reference beam indicating as expected 
that 0 T will be shifted in the same direction for both cases. Therefore the total 
shift in Of is the difference between the two: 


X|—X 2 k\— X 2 

a ,cos Oj a ,cos 0 r 

X|—X 2 fl|—fl 2 
cos Of a x ai 


From Bragg’s equation 

*1 

a,-:- 

2sin~{29—8) 

Xi 

or- 1 - 

2s/>iy(20-«) 

„ _ X,2»/i4(2W-6)-X 1 2sm^ 2 ®~«) 

o i Q 2 _ 1 2_ 2_ 
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~(2«uyU«)miy(S)]. 


-S*!*,,! 

Xj cos#/ 2 


Note that for a single reference beam, Le. 8-0, there is no shift in 0, for a change 
in X so double exposure single reference beam holograms are readily reconstructed 
using any laser wavelength without loss of fringe contrast. Again for the specific 
experiment to be discussed later, X|—532nm and Xj—514.5nm. For a fringe con¬ 
trast y-O.7, the uj is 5l4.5nm as before. Using Equation 26 it is found that S must 
be less than 3xl(T 3 degrees. Such a small separation between reference beams 
would cause overlap of the interfering image with the secondary images which 
were neglected in the development of Equation 17. These additional uncorrelated 
images add to the background illumination of the image field thus cutting the 
effective y in half and adding to the noise power of the detector. 
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Pulsed Holography 


Holographic recording using pulsed lasers has been used to reduce the effects 
of environmental vibration or object motion on holographic quality over a broad 
range of applications from holographic portraiture to nondestructive testing (Col¬ 
lier, 1971b). Nondestructive testing applications include contouring, vibrational 
mode analysis, and transient stress and deformation analysis (Gagosz, 1974; Chu, 
1972; Albe, 1976). In the area of transient analysis, pulsed holography has been 
employed to 'freeze' the motion of relatively large amplitude acoustic waves in¬ 
cluding plate waves and Rayleigh waves (Aprahamian, 1971; Gagosz, 1974). In 
both of these studies, however, the surface displacements were sufficient to be 
visible using single reference beam douhle exposure holographic interferometry. 
They are significant to this study in that they demonstrate that surface displace¬ 
ments, and rate of displacement in particular, associated with acoustic waves is 
insufficient to cause noticable reduction of holographic efficiency. Indeed the 
reduction of holographic efficiency may be computed directly for objects in uni¬ 
form and/or vibrational motion if the laser pulse duration, pulse shape, and the 
component of object velocity along the sensitivity vector is known (Gustafson, 
1980). Using Gustafson’s formulas, the reduction of holographic efficiency result¬ 
ing from the subfringe displacements encountered in this work would be less than 
0.5%. 
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Alt hough the reconstruction efficiency is not significantly reduced by small 
acoustic perturbations of the object surface, the accuracy with which such perturba¬ 
tions are recorded remains in question. Unfortunately the answer to this question 
is not available in the literature at present. However an argument can be present¬ 
ed based on the supposition that one cannot expect recorded phase accuracy any 
greater than that associated with the object displacement which takes place during 
the laser pulse. To gel a worst-case estimate of phase error using this argument it 
will be necessary to estimate the surface velocity normal to the object surface as an 
acoustic wave propagates across it. For the sake of simplicity, assume that the nor¬ 
mal (z) displacement, w(x,t) of the object surface resulting from a surface wave 
propagating in the x-direction may be described as a sinusoid as follows: 


where k a is the acoustic wave propagation constant 2 tt/A* i 
frequency. The displacement velocity, then is: 


-~^V0COs(aU—k a x), 


the maximum of which will occur when (a*— k a x W). Thus 


An acoustic frequency of 500kHz corresponds to an acoustic wavelength of about 
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6mm for aluminum where the surface wave velocity is 92% of the bulk shear wave 
velocity. In the experiments documented here, a frequency doubled NdiYAG 
laser (X—532nm) with a half amplitude pulse duration of 9ns is used. Therefore 
the phase uncertainty, A<£, associated with the finite recording pulse width is 

A«-2*(wow)r (27) 

- 2 S32 2 xV (2ff50 ° X 103)(9X l0 '* )wo 

—6.68x IOVq radians 
—3.83x 10 7 w 0 degrees. 


Thus for a peak wave amplitude wq— y or 266nm, the phase uncertainty A$—1°. 
Bear in mind that a surface displacement as large as 266nm will produce one holo¬ 
graphic fringe in the interference field. Consequently a 1° phase error would still 
permit measurement of the wave peak to one part in 360 or about 0.3%. Accord¬ 
ing to Equation 27, one may expect better base accuracy for smaller surface dis¬ 


placements. 
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Experimental Design 


Three primary experiments were performed to demonstrate that heterodyne 
analysis of pulsed holograms can be used to measure sub-fringe displacements as¬ 
sociated with traveling surface acoustic waves. In the first experiment, a CW 
heterodyne system designed very much after Dandliker’s was built. Using this sys¬ 
tem, the processing electronics were optimized and sub-fringe interpolation to 
1/900 of a fringe (2.9 a) was demonstrated. The last two experiments addressed 
the problems of synchronization of the laser output with the sound wave source 
for multiple exposure holography. As pan of the third and final experiment, holo¬ 
grams recorded using the frequency doubled NdiYAG laser were analyzd using 
heterodyne readout and phase detection. 


CW Heterodyne Holographic Interferometry 


The holographic recording system diagramed in Figure 14 was used to record a 
double exposure, dual reference beam hologram of the articulating surfaces of the 
tibia! component of an artificial lenee joint (see Figure 15). The two articulating 









'igure 14: CW recording set-up for 
heterodyne holography. 
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surfaces are concave and serve as sockets imo which the condyles of the upper, 
femoral, component are located. The object was placed in a sealed chamber such 
that the surfaces of interest could be viewed through the optically flat front win¬ 
dow (see Figure 16). A contour hologram of the object was recorded using the 
multiple refractive index technique referenced earlier. By changing the pressure 
within the test chamber between exposures, the contour interval could be continu¬ 
ously adjusted so that less then one contour fringe could be observed over the im¬ 
age. Thus the details of the surface contour were revealed only by sub-fringe inter¬ 
polation. 

When recording a hologram, the beam from a Spectra-Physics 165 argon ion 
laser was split and directed through three acousto-optic modulators (refer to Figure 
14). Each modulator (Isomet 1205C) was driven by an 80MHz signal from a sin¬ 
gle oscillator/driver during exposure. Beam ratios were adjusted so that either 
reference beam was one to two times as bright as the object beam at the film 
plane. A different reference beam was blocked for each exposure. Agfa 10E56- 
N AH film in a 4"x5" glass plate format was used. Total exposure was approximate¬ 
ly 10 ergs/cm 2 . The film was developed for 5 minutes at room temperature using 
Kodak D-19 developer. After a 30 second rinse, the plate was fixed with Kodak 
rapid fixer for 5 minutes. 

Once dry, the film plate was replaced in its holder for reconstruction and 
analysis. The beamsplitter used to provide illumination of the object was removed 
as was the etalon within the laser cavity. Both these steps increased the amount of 
light available in the reference beams. Precise replacement of the film plate was ac- 









Figure 15: Artificial knee implant. 











Figure 16: Test chamber for contouring 
knee implant. 
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coraplished by observing the interference pattern in the film plane and adjusting 
the film position and final reference mirrors to minimize the fringe order. Once the 
film plate was repositioned, the drive signal from one of the acousto-optic modula¬ 
tors was changed from the common 80MHz oscillator to a second oscillator this 
time running at 80.1MHz. Merit model R1500 silicon PTN diodes with integral 
transconductance amplifiers were used for signal detection. The reference detector 
was positioned for best S/N ratio. A second detector was mounted on an x-y 
translation stage driven by Newport Research model 850-1 linear actuators. The 
reference and probe signals were amplified by Tektronix AM501 differential 
amplifiers, passively hand-pass filtered at 100kHz and fed to the reference and sig¬ 
nal inputs of a PAR model 5301 lock-in amplifier. The lock-in amp was configured 
to measure phase difference between the reference and probe signals. The analog 
phase output signal was connected to the Y-axis input of a Hewlett-Packard 7045 A 
x-y recorder. The x-axis was driven with a ramp signal from an HP 3304 function 
generator. (In subsequent experiments probe detector positioning, phase voltage 
recording, processing, and display were all performed under computer control. The 
implementation of this computerized data acquisition and analysis system will be 
discussed later.) Using the actuator joystick controller, the probe detector was po¬ 
sitioned to one side of the output image and then scanned across the image at a 
constant velocity. The beginning of the probe scan was synchronized with the 
recorders x-axis drive signal. Ry incrementing both the recorder pen and the 
detector probe along the y-axis, a y-moduiated raster scan of the knee prosthesis 


was generated. 
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Double Pulsed, Single Reference Beam Holography 

A Quantel model 580-C NcLYAi} Laser with an etalon and frequency doubling 
crystal was used to record all of the pulsed holograms. The laser has a horizontally 
polarized output with a pulse width if 9ns and an output energy of 1.2 joules/pulse 
at 1064nm and about half that en> rgy at 532nm. The laser energy in the green 
(53*-tm) can be reduced either by lentional misalignment of the doubling crystal 
or by disabling the amplifier flash .tips. In the later case, the maximum energy 
per pulse at 532nm from the osc itor alone is 350mJ. It was this low evergy 
configuration which was used to n ird holograms of specularly reflecting objects. 
Recording of diffuse object surface .squired higher illumination intensities. 

The recording geometry diagrar d in Figure 17 was used to record the travel¬ 
ing surface wavefront propagating ac >ss an aluminum plate. The plate dimensions 
were 25cm x 38cm x 3.8crn thick. Tne front surface of the plate was grit blasted 
to provide a diffusely reflecting surface. The glass plate beamsplitter was posi¬ 
tioned near the brewster angle to permit most of the beam energy to pass through 
it to illuminate the object. 

The purpose of this experiment was in part to verify the timing and allignment 
of the pulsed holographic recording system. Since no fringe interpolation was to 
be done at this stage it was necessary that the recorded acoustic wavefronts be 








-75- 







- 76 - 


sufficienily large to produce interference fringes which may be visually observed 
when the hologram is reconstructed- In order for this to occur the the surface wave 
displacement amplitude must be greater than A/2 or about 2500A. To generate 
such large surface displacements a small explosive charge was used as the sound 
source. The explosive used was silver acetylide which is a white precipitate formed 
when acetylene gas is bubbled through an aqueous solution of silver nitrate. Once 
rinsed and dried, the silver acetylide powder was pressed into small pellets about 
3mm in diameter and 0.5mm thick. Each pellet was then sandwiched between 
layers of cellophane tape and taped to the surface of the plate at the leftmost edge. 
The infrared output of the pulsed laser was focused on the pellet to detonate the 
explosion. 

Assuming that the surface wave velocity is about 92% of the shear wave veloci¬ 
ty in aluminum (3100 m/s), then the surface wavefront takes about 45 ms to pro¬ 
pagate half way across the plate surface. Since the same laser was used to detonate 
the explosion and record the hologram, it was necessary to modify the pockels cell 
triggering circuit to double pulse the laser during a single flash lamp pulse with a 
45ms pulse separation. The pockels cell trigger circuit (Figure 18) uses eleven 
transistors in series which operate in avalanche mode to deliver a high voltage 
from capacitor C| to the pockels cell. Once avalanche is initiated the voltage 
across the series transistors must fall well below the avalanche threshold voltage 
before they will again turn off. It is the time required for this voltage decay which 
limits the rate at which the pockels cell may be triggered. To decrease the time re¬ 
quired to reset the series transistors, an additional transistor was inserted across 
the pockels cell output terminals. This transistor turns on shortly after the high 
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Figure 18 


: Q-switch trigger circuit. 
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voltage is applied to the cell Once on, it provides a low resistance path which 
more rapidly bleeds charge from the capacitor thus causing a correspond^ in¬ 
crease in the voltage decay rate. Once implemented, this modification permits 
double pulsing of the laser with pulse separations externally controlled over a 
range from 20/xs to 100ms. 

Optimum double pulsing of the laser required proper timing of the Q-switch 
pulses relative to the beginning of the flashlamp pulse. For this reason, the trigger 
system shown in Figure 19 was used. An adjustable digital delay generator (BNC 
model 7010) was used to trigger an HP model 8005R pulse generator with double 
pulse output. The delay unit was triggered by the flashlamp drive signal generated 
internally by the laser control electronics. A simple buffer circuit capable of provid¬ 
ing the required 15 volt Q-switch trigger pulses into the 50 ohm trigger input cir¬ 
cuit was designed and built since the TTL output of the HP pulse generator was 
insufficient to trigger the Q-switch. 

The sequence for recording the pulsed acoustic surface wavefront was relatively 
straight forward. The pulse separation was set for 45 ms and the digital delay adjust¬ 
ed to provide equal intensities in both pulses. The explosive charge was fixed to 
the plate surface and a film plate placed in the holder. When the laser was fired, a 
holographic exposure of the plate surface was recorded using the green line of the 
laser output. At the same time, the light from the infrared port of the laser de¬ 
tonated the explosive charge. Forty-five microseconds later, the laser emitted a 
second pulse of light. This time a second exposure of the hologram was made us¬ 
ing the green wavelength so that any surface deformations caused by the explosion 








Q-SWITCH TRIGGER 
FROM LASER 



TO Q-SWITCH 
TRIGGER CIRCUIT 



Figure 19: Laser trigger system 
and timing. 
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were recorded. The infrared beam again struck the plate. However, it was of no 
consequence since this time there was no charge to be detonated. Again Agfa 
10E56 film was used for recording. Kodak D-19 developer for 5 minutes was fol¬ 
lowed by a water rinse. Instead of fixing these plates, however, the plates were 
beached for about 2 minutes in a solution of potassium dichromate and sulfuric 
add to improve diffraction efficiency. 


Surface Acoustic Wavefront Recording for Heterodyne Analysis 


Through this experiment, the attempt to measure and display subfringe detail 
of displacements associated with surface wave propagation was made. The optical 
arrangement diagrammed in Figure 20 was used to record dual reference holo¬ 
grams of the surface wavefront generated by detonation of a silver acetylide pellet 
as described in the previous experiment. The same aluminum spedmen used in 
the previous experiment was again used as the test object. This time however, the 
object surface was polished to reduce scattering of the collimated illuminating 
beam. 

To fadlitate heterodyne analysis, a dual reference beam double exposure holo¬ 
gram could have been recorded as described previously in the Background section. 
Using a double-pulsed laser to both detonate the expolsive charge and expose the 







DIVERTING PRISM 


Figure 20: Triple exposure recording 
geometry. 
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holograms, however, would have required some type of electro-optic switch to 
change the reference beam path during the 45*ts delay between the laser pulses. 
Unfortunately, such a switching device was not available for this experiment. In¬ 
stead, a dual reference beam, triple exposure hologram was recorded. The first ex¬ 
posure was made using reference path 1 as shown in Figure 20. During this first 
exposure, the IR port of the laser was closed to prevent detonation of the explo¬ 
sive charge. After the diverting prism was removed, the IR port was again opened 
and the laser was double-pulsed to record exposures 2 and 3 using reference path 
2 . 


While the triple exposure technique permits heterodyne analysis of the holo¬ 
graphic interferogram, it has a pronounced effect on the calibration and sensitivity 
of the analysis. This effect can be understood by examining the mathematical 
description of the output of the tr ole exposure dual reference hologram. The in¬ 
tensity distribution in the output p one of the interferometer is nearly identical to 
that described by Equation 17 exce pt that an additional image field, Oj, must also 
be added in. Thus the intensity in the output plane becomes 

/= | o ]e J ' m " ,+ *' (x ' y)) + 0 p J(m l , ' h *{ x ' y)) +i ) jg jU .•+4jt*0'W |2 

■"of-t-of +20j03C0s(«ii 1 —0 2 >-(-2OiO2COS(fc> 1 —«U 2 )r-f20i03C0s[(«i— a^Ir-K^l - ^)] 

—o j 1 +2o f +2o f cosA<£-K4o iO 2 COS-^p )cos(^ar ^-—-) 


for 
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OjP-03 

A <^0|—<^2 

Aai—<i>2 

/-Afcos(Aa><4~^-) 


where 


. A <4 

4o,02cos-r- 

_ 2 

r» 2 4-2o f (1+coa A«£) 

It is clear from the final form of this expression that the measured phase 
differences using, the triple exposure technique will be half that which would be 
measured using the double exposure technique. Further, the modulation index 
will be zero when A <4 has a value equal to odd multiples of 180° and will be nega¬ 
tive for 180°<A<4<540° the effect of this negative modulation index is to generate 
a 180° discontinuity in the output of the phase detector when the optical phase 
difference between corresponding points on 0 \ and o 2 exceeds 180°. Fortunately, 
for subfringe displacements, optical phase differences will always lie between 
4-180° and -180° thus insuring that the modulation index will always be positive. 
As the modulation index falls below 1, the signal to noise ratio degrades as 
descnbed earlier. Again for small displacements, this may not be a real problem 
since M will always be between 0.8 and 1.0 for optical phase differences between 
0° and ±120°. 


Talcing the above into consideration it was felt that the triple exposure tech¬ 
nique could be used to demonstrate the application of heterodyne analysis to 











pulsed holographic interferometry. The triple exposure hologram was recorded on 
Agfa 10E56NAH plates and processed as before The processed plates were then 
placed in the heterodyne readout system diagrammed in Figure 21. The plates 
were reversed in the film holder so that the conjugate, real image was projected on 
the output plane without requiring imaging lenses. Alignment was performed with 
both A-O modulators driven at 80.0MHz so that a stationary fringe pattern could 
be observed. Once the system was aligned, the modulators were driven at 
80.0MHz and 80.1MHz, respectively. 

Figure 22 shows the data collection system in block diagram form. Again, two 
Merit R1500 optical detectors were used to pick up the 100kHz beat signal. This 
time, however, no beamsplitter was used to provide a reference signal. Instead, 
the reference detector was placed in the output plane at one corner of the project¬ 
ed image. The scanning detector was mounted on an x-y translation stage driven 
by two NRC 85C-C motorized micrometers. The micrometers were driven through 
an interface controlled by a Rockwell AIM 65 microcomputer. AH micrometer 
drive routines were written using FORTH programming language and are listed in 
the Appendix. After selecting the appropriate vertical position with the y-axis mi¬ 
crometer, the horizontal scan routine was run. This routine scanned the photo- 
detector over a distance of about 1cm in 125 fim increments. After each increment, 
a trigger pulse was sent by the AIM 65 to the A/D converter module of an IBM 
model 7406-1 Device Coupler. The trigger pulse signaled the module to convert 
the analog phase difference signal from the lock-in amp and to transmit that digi¬ 
tized value through an RS232 interface to an IBM 5100 computer. The 5100 com¬ 
puter was used to decode the digitized data and store the data in vector form such 
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that each 90 element vector corresponded to one scan across the output field The 
raw data received by the computer is modulo 360° in that phase values measured 
by the lock-in amp are restricted to those between -180° and +180°. Consequent¬ 
ly, a computer routine was used to "unwrap" the raw data to remove these large 
discontinuities. As predicted with the relatively large surface displacements associ¬ 
ated with an explosive sound source, 180° discontinuities in the data are also ob¬ 
served These too are removed by computer processing. Finally the processed 
data was plotted on a Tektronix 4013 graphics display terminal and subsequently 
plotted on a Hewlett-Packard 7045 A x-y recorder using the Analog Output Module 
of the Device Coupler. All collection and processing software in the IBM 5100 
computer was written in APL and is listed in the Appendix. 






Results and Discussion 


Sub-fringe Interpolation: Contour Holography 


The image shown in Figure 23 is the output of the dual reference contour 
hologram of the knee implant recorded as described in the previous chapter. In 
order to photographically record the contour fringes in the output field, both A-O 
cells in the reference legs were driven at the same 80.0MHz drive frequency. Dur¬ 
ing heterodyne analysis, the two modulating frequencies differ by 100kHz thus 
causing the fringe pattern to fluctuate at such a high rate that the fringes are no 
longer visible. Optical detectors in the output plane, however, detect the 
sinusoidally varying intensity and transduce it into a voltage signal analog. Initial¬ 
ly, the position of the scanning detector relative to the reference detector in the 
output plane was fixed and the phase difference as a function of time was plotted 
in order to measure system noise. Using this technique, the short term phase 
noise was characterized as having peak-to-peak fluctuations of 0.288°. Long term 
drift of about 0.04° per minute was also observed. In subsequent experiments, the 
drift rate was observed to deviate from this value indicating that drift may be in 
pari the result of dimensional changes in the apparatus in response to changes in 
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environmentai temperature. 

The entrance pupil of the imaging optics was the film plate itself. Therefore 
the speckle size was determined by the 4 inch (10cm) dimension of the film plate. 
Given that the image plane was some 50cm behind the film plate, the speckle di¬ 
ameter was about 2.6 microns. The number of speckles, N, within the 1mm diam¬ 
eter detector was then about 1.5x10 s . From Equation 23, the expected phase error 
with iV— 1.5X10 -4 and an estimated modulation index of 0.8 was about 0.11°. 
Therefore the overall phase noise may be as large as the sum of the elctronic and 
speckle noise or 0.4°. It should be noted that this 0.4" noise level corresponds to 
an optical path length sensitivity (for SNR — 1) of 0.57nm (5.7A) or an object 
displacement as small as 0.29nm (2.9A). In the specific case of contour hologra¬ 
phy, the sensitivity is 1/900 of the contour fringe interval. 

A single horizontal line scan was made scanning from left to right across the 
image plane and passing through what would be the center of the stationary fringe 
pattern shown in Figure 23. Note that the plot of relative phase difference as a 
frunction of position is modulo 360" (Figure 24) and that if "unwrapped" it would 
provide a profile of the curved surface. However, sufficient sensitivity existed us¬ 
ing this analysis system to detect the surface features of interest e'ren when the 
holographic contour interval was larger than the depth of the surface relief over 
the entire object. Therefore, a second contour hologram was prepared using the 
double refractive index technique where, as described previously, only a slight 
change in air pressure inside the lest chamber was made between exposures. The 
resulting contour image showed less than one complete fringe over the object. 
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Figure 24: Single line scan across knee implant. 
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Thus when ynnit^ over the image, the phase difference between the s ca nn i n g 
arri refereiKc signals never exceeded 360°. Figure 23a shows a series of horizontal 
over a section of the image. The spacing between the scans corresponds to a 
vertical displacement of about 0.5mm along the image. The scans follow the nor¬ 
mal curvature of the artificial knee socket. Just short of hallway across, all but the 
lowest scans show a deviation from the smooth contour. This deviation was 
caused by the gouge placed in the socket which was not apparent in the conven¬ 
tional contour hologram (Figure 23). A section of the gouge is shown in greater 
detail in Figure 25b by icreasing the gain of the chart recorder. The displacement 
between scans is 0.1mm in this figure. 


Double Pulsed, Single Reference Beam Holography 


Photos of thcj reconstructed images of the aluminum plate recorded as 
described in the corresponding section of the previous chapter are shown in Fig¬ 
ures 26 a, b, and c. Figure 26a shows the acoustic wavefront as it appeared 35 mi¬ 
croseconds after the explosive charge was detonated. Figures 26b and 26c similarly 
show wavefronts from subsequent experiments where the laser pulse separation 
was 40 and 45 microseconds respectively. 
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Frora these photographs, the wave velocity and therfore the wave type can be 
determined. Table I summarizes the velocity data obtained from Figures 26 thus 
confirming that the observed primary wavefronts travel at a velocity corresponding 
to the Rayleigh wave velocity. Following the primary wavefront, one observes an 
irregular pattern of fringes. This pattern varies from experiment to experiment as 
the laser pulse interval is varied. Presumably this pattern is the result of surface 
displacements caused by the superposition of waves which are multiply reflected 
from the front and back surfaces of the test plate. 


Heterodyne Analysis of Surface Acoustic Wavefronts 


From the previous experiment it was confirmed that the primary surface wave- 
front passes across the center of the test plate some 45 microseconds after the ex¬ 
plosive charge is detonated. This single data point was important since those holo¬ 
grams recorded for subsequent heterodyne analysis were made using the polished 
surface of the specimen as discussed earlier and only an 8cm diameter region was 
illuminated by the object beam. It was important, therefore, that the wavefront be 
passing through the 8cm circle when the final exposure of the hologram was made. 


Two steps were involved in the alignment of the developed holograms in the 
heterodyne armlysis system. Driving both A-0 cells at the same frequency, a video 
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caraera was focussed on the image plane. A cross which marked the center of the 
test surface was visible in the images reconstructed by both reference beams. 
Coarse alignment of the two images was obtained by changing the two reference 
beam angles to superimpose the two crosses observed on the video monitor. Final 
alig n ment was made by focussing the camera on the plane of the hologram plate 
and adjusting the reference beam angles to minimize the order of the observed 
fringe pattern. As predicted, the 10° separation of the reconstructing beams and 
the use of different recording and playback wavelengths made it impossible to el¬ 
iminate all fringes observed in the hologram plane. Under the best alignment con¬ 
ditions, less than one fringe was observed over an area about 3cm in diameter in 
the center of the illuminated region. Beyond this central region, multiple fringes 
encircling the center were observed. Focussing back on the image plane, the effect 
of wavefront misalignment manifested itself in reduced fringe contrast and fringe 
distortion beyond the 3cm central region. Photographs of the output plane with 
proper and improper beam alignment are shown in Figure 27. Note the distortion 
of the fringes in Figure 27b resulting from slight misalignment. 

Figure 28a shows relative phase difference as a function of position across a 
lcm length of the image. This raw data is once again modulo 360* and exhibits 
the 180° discontinuities predicted for triple exposure heterodyne holography. Ac¬ 
counting for the discontinuities and "unwrapping" the data produces the plot of dis¬ 
placement versus position of Figure 28b. Two relatively narrow peaks are ob¬ 
served in the leading portion of this wavefront plot. Although it was beyond the 
scope of this research to explain these peaks, it was originally thought that they 
may be the result of the superposition of the surface wave emanating from the 






















a) Properly aligned 



b) Inproper alignment 


Figure 27: Image plane recordings of dual reference 
beam holograms. Note curvature of fringes 
in improperIt aligned case. 
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source location with the wavefront which is reflected from the edge of the test 
plate immediately to the left of the charge. To test this hypothesis, a second ex¬ 
periment was performed identical with the first except that the charge location was 
displaced towards the center of the piate a full centimeter from the plate edge. As¬ 
suming that the twin peaks were the result of interference between the direct and 
reflected surface waves, the wave shape and especially the separation between the 
peaks should have changed. Instead, as can be seen from Figure 29, both the 
peaks and their relative spacing remained the same even after the charge location 
was changed. Therefore, edge reflection does not account for the dual peaks ob¬ 
served in the surface wave shapes. By changing vertical positions between scans, a 
full field raster over an area of the test surface can be obtained and plotted in *3- 
DT as shown in Figure 30. 

Several sources of inaccuracy limited the resolution of this heterodyne analysis 
system. First, the measured efficiency of the triple exposure holograms was less 
than 1% providing an output field intensity of about 1 microwatt/cm 2 . From 
Equation 12, the expected phase variation was 0.06" for a fringe contrast of 0.7, a 
lms measurement interval, and a 10kHz bandwidth. In fact, the measured signal 
variation for a stationary scanning detector was 0.6 s , which was much worse than 
predicted. The low signal intensity in the output plane and the relatively poor 
noise figure of the amplifiers used would tend to indicate that the detection system 
was not shot noise limited (a key assumption in deriving Equation 12). In addition 
some of the noise may be attributed to variation in pathlength owing to environ¬ 
mental effects. 
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Figure 29: Displacement plot fox 
image recorded with t 
by 1cm towards the ce 



















A second noise source resulted from the effort to conserve light in the output 
plane by eliminating the need for focussing lenses. As a result, the distance, z, 
between the hologram and the output plane was a fairly long 61cm. The effective 
aperture of the imaging system was the 8cm diameter exposed region of the holo¬ 
gram. The speckle size was determined to be about 4.8 microns in diameter so 
that 4.4x10* speckles were encompassed by the 1mm diameter detector. Using 
Equation 23 the phase error resulting from speckle noise was on the order of 
0.28°. 

Finally, the rather large diameter of the detector itself caused some spatial 
averaging or smoothing of the surface features. As a result, one can not be certain 
of the absolute height of the leading edge peaks observed in Figures 31 , 32 , and 
33. It is quite possible that these peaks are considerably larger than indicated. 




Conclusions 


Subfringe detail of traveling surface acoustic waves may be measured and 
displayed in full field using heterodyne holographic analysis of multiple exposure 
holograms recorded with a pulsed laser. Systems for recording and displaying sur¬ 
face wavefronts have been demonstrated to provide displacement resolution to 
0.43nm (4.3A) over a 3cm diameter field. The limited field size is the result of 
using different recording and playback wavelengths. The field size may be im¬ 
proved to the size of the entire field recorded by the hologram by changing the 
recording and playback geometries to make the two reference beams more nearly 
colinear. Indeed when the same wavelength was used for both recording and 
analysis of the contour hologram, analysis of the entire field was possible. 

The single most important step m improving the sensitivity of the heterodyne 
analysis system would be to increase the intensity of light in the output field. 
Several options are available to increase output intensity. First, holographic 
efficiency may be greatly improved by using a different recording medium either 
during the recording of the hologram itself or by copying the primary hologram 
onto a second medium using a longer CW copying process. Dichroma ted gelatins 
and photoresist holograms have been reported to have efficiencies over 90%. Using 
an electro-optic switch, as described earlier, to allow double rather than triple ex¬ 
posure holograms to be recorded for heterodyne analysis would not only improve 
holographic efficiency but would provide a constant fringe visibility in the output 
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plane independent of optical phase differences. Still more improvement could be 
expected if a higher frequency detector capable of sensitive recording at up to 
80MHz were used. This would permit the use of only a single A-0 cell for fre¬ 
quency shifting. Using the undiffracted beam for one reference, and the diffracted 
beam for the other reference, virtually all of the light from the laser would be pro¬ 
jected on the hologram. Presently, the undiffracted beams from the two A-0 
modulators are blocked off and that energy is lost. 


Increasing the light available in the output plane would also make it possible to 
further reduce spatial noise in the system. With more light available, a lens system 
may be used to reduce the working distance behind the hologram and thereby 
further decrease speckle size and thus improve phase accuracy. In addition, the 
detector area may be reduced to further improve lateral resolution as well. 


The results of this work appear to support hypotheses proposed by the author 
and others working in the field of heterodyne holographic interferometry. Rased 
on the experimental results obtained and the mathematical relationships which 
have been developed to describe system performance, the author feels confident in 
predicting that a practical system for studying sub-Angstrom transient phenomena 
may be realized using heterodyne holographic interferometry. 
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